arXiv: 1502.01189vl [gr-qc] 4 Feb 2015 


The universe dominated by the extended Chaplygin 
gas 


E.O. Kahya,® B. Pourhassan 6 

“ Physics Department, Istanbul Technical University, Istanbul, Turkey 
b School of Physics, Damghan University, Damghan, Iran 

E-mail: eokahya@itu.edu.tr, bpourhassan@yahoo.com 

Abstract: In this paper, we consider a universe dominated by the extended Chaplygin 
gas which recently proposed as the last version of Chaplygin gas models. Here, we only 
consider the second order term which recovers quadratic barotropic fluid equation of state. 
The density perturbations analyzed in both relativistic and Newtonian regimes and show 
that the model is stable without any phase transition and critical point. We confirmed 
stability of the model using thermodynamics point of view. 
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1 Introduction 

Accelerated expansion of universe can be described by dark energy [1, 2], In that case there 
are several dark energy candidates. One of the first and simplest model is the cosmological 
constant [3, 4] which has two important problems known as fine-tuning and coincidence. 
Therefore, alternative models were investigated such as quintessence [5-11] and k-essence 
[12, 13]. 

Another attempt to solve the coincidence problem is interacting dark energy model [14- 
23] . The dark energy density is of the same order as the dark matter energy density in the 
present universe, therefore one could imagine that there is some connection between dark 
energy and dark matter. 

The above models are based on modification of the scalar field and the fine-tuning of the 
potential. Instead, one can modify the equation of state to produce an exotic background 
fluid so-called Chaplygin gas. The pure Chaplygin gas emerged for the first time from the 
string theory point of view [24, 25], then it is used to study classical solution of bosonic d- 
brane in {d+ 1,1) space-time [26]. One can begin from Nambu-Goto action in the light-cone 
gauge and solve the momentum constraint to obtain the Chaplygin gas equation. Then it 
was seen that different kinds of matter on the brane, e.g. Chaplygin gas, can be introduced 
[27]. The Chaplygin gas equation of state is given by [28], 

B 

PCG = -, (1-1) 

PCG 

where B is an arbitrary parameter which usually considered as a constant. In that case 
density perturbations in a universe dominated by the Chaplygin gas has been studied by 
the Ref. [29]. 
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In order for the universe to pass from a dust dominated epoch to a de Sitter phase through 
an intermediate phase described as mixture of cosmological constant and radiation the 
Chaplygin gas equation of state generalized to [30, 31], 

B . . 

PGCG = —s-’ ( L2 ) 

Pgcg 

where a and B are free parameters. It is clear that a = 1 reproduces the pure Chap¬ 
lygin gas model. This model is called the generalized Chaplygin gas model [32]. In the 
recent years the generalized Chaplygin gas was the subject of several cosmological and 
phenomenological studies. For example Brane-world cosmology and Brane-world inflation 
including the generalized Chaplygin gas has been studied [33, 34], There is also possibility 
to consider viscosity in the generalized Chaplygin gas model [35-39]. 

There is also a class of equations of state that interpolate between standard fluids at 
high energy densities and generalized Chaplygin gas fluids at low energy densities which 
is called the modified Chaplygin gas [40]. Viscous modified Chaplygin gas was introduced 
by the Ref. [41] and completed by the recent works [42-44]. Modified Chaplygin gas was 
introduced with the following equation of state [45, 46], 

. B . . 

PMCG = ApMCG - ^ -’ (l- 3 ) 

Pmcg 

where A, a, and B are parameters of the model. The case of A = 0 recovers generalized 
Chaplygin gas equation of state, and A = 0 together a = 1 recovers the original Chaplygin 
gas equation of state. The MCG equation of state has two parts, the first term gives an 
ordinary fluid obeying a linear barotropic equation of state, while there are some models 
with quadratic equation of state [47-49], 

p = p 0 + u>ip + uj 2 p 2 , (1-4) 

where po, u\ and ui 2 are constants. Easily we can set po = u 2 = 0 to recover linear 
barotorpic equation of state. 

Modified Chaplygin gas include only barotropic fluid with linear equation of state, while 
it is possible to extend them to including quadratic barotropic equation of state given by 
(1.4). It yields to introducing extended Chaplygin gas with the following equation of state 
[50-53], 

PECG = ^ AnPECG a • (l- 3 ) 

Pecg 

It is obvious that the n = 1 reduced to the modified Chaplygin gas. In this paper, we would 
focus on the second order term which recovers quadratic barotropic equation of state, 

PECG = AipECG + MPECG - a -’ (l- 6 ) 

Pecg 

where a, A \, A 2 and B are free parameters of the model. 

On the other hand the temperature behavior and the thermodynamic stability of the gener¬ 
alized Chaplygin gas has been studied by the Ref. [54] , and it is found that the generalized 
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Chaplygin gas cools down through the expansion without facing any critical point or phase 
transition. Also, thermodynamics of the generalized Chaplygin gas was investigated by in¬ 
troducing the integrability condition, and thermodynamicalquantities have been derived as 


functions of either the volume or the temperature [55]. Validity of the generalized second 


law of gravitational thermodynamics in a non-flat Friedmann-Robertson-Walker (FRW) 
universe and an expanding Godel-type universe containing the generalized Chaplygin gas 
was confirmed by the Refs. [56] and [57] respectively. 

In extension of the Ref. [54], the similar work was performed for the case of the modified 
Chaplygin gas [58] and the same result was obtained. More discussion on thermodynam¬ 
ical behavior of the modified Chaplygin gas can be found in the Ref. [59]. Also, Ref. 
[60] extended the work of [55] to the case of the modified Chaplygin gas. Validity of the 
generalized second law of thermodynamics in the presence of the modified Chaplygin gas 
was investigated by the Ref. [61] and it was observed that the generalized second law of 
thermodynamics always was satisfied for the modified Chaplygin gas model. The gener¬ 
alized second law of thermodynamics in the brane-world scenario including the modified 
Chaplygin gas was verified for the late time behavior of the apparent horizon [62], 

In this paper, we study a cosmological model where the universe is dominated by the 
extended Chaplygin gas. Indeed, we would like to investigate the evolution of density per¬ 
turbations, trying to verify if the dominance of the extended Chaplygin gas is compatible 
with the formation of the universe. In that case we consider both Newtonian and relativis¬ 
tic descriptions of the extended Chaplygin gas. We find that the universe dominated by 
the extended Chaplygin gas is stable without any phase transition and critical point. We 
also investigate the thermodynamical aspects of the model and confirm the stability of the 
model using laws of thermodynamics. 

2 Extended Chaplygin gas cosmology 

Einsteins equations are given by, 




( 2 . 1 ) 


where the energy momentum tensor is 


T^u = {p + p)u^u v - pg^u- 


( 2 . 2 ) 


Assuming flat FRW universe, 


ds 2 = dt 2 — a 2 (dr 2 + r 2 (dd 2 + sin 2 ddc/) 2 )), 


(2.3) 


yield to the following Friedmann equations [63] 



H = ~~{PECG + PECG), 


(2.4) 
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where H is the Hubble expansion parameter and a is the scale factor. Then conservation 
equation read as, 

Pecg + 3-{pecg + PECO ) = 0, (2.5) 

a 

where 8-/r G = 1 is used. In order to simplify the calculations and to reduce the number 
of free parameters of the model and also to obtain an analytical expression of the energy 
density in terms of the scale factor, we assume the following conditions, 

a = 1 , 

A\ = A 2 — 1, 

B = 2 A 2 . (2.6) 

So, the equation of state (1.6) yields to the pure Chaplygin gas plus quadratic barotropic 
fluid equation of state. While it is very special case of the extended Chaplygin gas, we will 
show that the model may be suitable to describe universe. Having an explicit expression 
of the energy density in terms of the cosmological parameter is our motivation to choose 
the above condition. Therefore, only free parameter of the model is A 2 , and we can solve 
the conservation equation (2.5) using (1.6) and (2.6) to obtain the following relation, 

, = ln (pecg + 2 PECG + 2) _ In ( pecg - 1) _ arctan ( p E CG + 1) „ (2?) 

30 A 2 15 A 2 5 A 2 ’ [ ' ’ 

where C is an integration constant. In order to obtain an analytical solution we use 
tan -1 ( pecg + 1) « 7r/2 approximation which is exact for for pecg ^ 1- In the Fig. 
1 we show that exact and approximate solutions have similar manner and we can use 
approximate solution, so we can rewrite the equation (2.7) as follows, 



Figure 1 . Typical behavior of p in terms of a. blue: exact solution, red: approximate solution. 


i 2 +V-l + 5 a 30A2 e 37r 

Pecg =1+ 


1 + a 30A2 e 3ir e 37T ’ 


( 2 . 8 ) 









where tan -1 ( Pecg + 1) ~ 7t/ 2 approximation is used. This approximation is valid when 
Pecg 1 corresponding to the early universe. However our solution will be valid at all 
times and our approximate solution is very close to the late time behavior with pecg ^ 
1. This is due to the fact that tan^ 1 ( pecg + 1) ~ 7t/ 4, for pecg 1 therefore small 
compared to the logarithm term. 

Hence, the Hubble parameter can written as follow, 


H = 



2 + V-l + 5a 30A2 e 37r 
-3 + 3a 3(M2 e 37r 


(2.9) 


Solving the above equation for a gives the scale factor in terms of time. We can obtain the 
following relation, 

f — /(o) + g(a) + C = 0, (2.10) 

where C is an arbitrary integration constant, and we defined, 


/(«) 

s(«) 


V31n (2 V-l + 5e 37r a 30A2 + l) 
15A2 

2\/3arctan (V —1 + 5e 37r a 30A2 ^ 

15A 2 


( 2 . 11 ) 


The contribution of the second function g(a) is small compared to the logarithm term, so 
we can remove it and fix the constant C to capture the effect of g(a). Therefore, we can 
obtain the following time-dependence for the scale factor, 


g 10 a 2 


a = 


20 


[X{tf - 2X(t) + 5] 30A 2, 


( 2 . 12 ) 


where we defined, 


X(t) = e 5 ^ A ^ t+t °\ (2.13) 

and g{a) + C = to is used. It is clear that the scale factor is an increasing function of time. 


3 Relativistic analysis 

We can study density perturbations of the Einsteins equations using the following relations, 

9fiu = 9nv T 

p = p + 5p, 
p = p + 5p, 

u,j. = u^ + Sufj,, (3.1) 

where h^ v ,5p,5p, and 5u^ are small perturbations around g lw , p,p-, and respectively. 
Assuming plane wave expansion to describe the spatial behavior of the perturbations and 
also the synchronous gauge condition h^o = 0 give us the following general equations [29], 

h + 2Hh = 5p + 35p, (3-2) 
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(3.3) 


5p+(p + p)6 + 3(p + p)5p ~^(p + p)h = 0, 

71 ? 

(p + p )0 + (p + p)0 + 5 H(p + p)9 = —^5p, 

a z 


(3.4) 


where n is the wavelength of the perturbations, 6 is velocity potential with definition 
V.;0 = Ui [26] and 5 = 5p/p. Also, 



(3.5) 


We can obtain behavior of 5 using simple assumption that 6 be a constant. In that case 
from the equations (3.3) and (3.4) we can obtain, 


. T 7 (p5 + Sp - 2H5p) - 2 c(p + p) - 5 cH(p + p) 

6 = - -;- 


(3.6) 


p + p + 5H (p + p) + 5H(p + p) 



(3.7) 



+ (P + p)c +(p + p)0 - ?—^-\p6 + 3 p{p + p)6 + (p + p)9], 


(3.8) 


p + p 


where 9 = c (c is a constant) is used. Now, using the equations (3.6), (3.7), (3.8) in 


the equation (3.2) give us a second order differential equation for 5. Note that using the 
equation (2.12) we can obtain p, p and H in terms of time. Resulting equation can be 


solved numerically to give us the behavior of 5 as illustrated in the Fig. 2. We can see that 
perturbation may vanish at the late time for n ~ 3. In order to compare the result with 


the previous versions such as modified Chaplygin gas we assume an infinitesimal value for 
A 2 and obtain results which are presented in the Fig. 3. We can see opposite behavior, 
so perturbations may vanish for n = 0 and initially it behaves as t m (m = 2/3 for dust), 


which agrees with the result of the Ref. [29]. As we know, in the limit of small wavelength 


of the perturbations, the relativistic problem reduces to the Newtonian problem. So, it is 
interesting to study the perturbations in Newtonian description which we will perform in 
next section. 

4 Newtonian description 

The evolution equations for perturbation in Newtonian description are given by [29], 



where, 



(4.2) 


is sound speed which is always positive for A 2 > 1 or A 2 < 0. Later we will give more 
details about the squared sound speed. In Fig. 4 we can see that perturbations vanish at 
the late time and the model is completely stable without any phase transition and critical 
point. 
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Figure 2. Typical behavior of S in terms of t for finite Ai (extended Chaplygin gas), n = 0 (blue), 
n = 1 (green), n = 2 (yellow), n = 3 (red). 



Figure 3. Typical behavior of S in terms of t for infinitesimal A 2 (modified Chaplygin gas), n = 0 
(blue), n = 1 (green), n = 2 (yellow), n = 3 (red). 


5 Thermodynamics 


The beginning point of thermodynamical behavior of the model is the relation of some 
thermodynamical quantities with the energy density [60], 


P = 


U_ 

V’ 


(5.1) 


where U is the internal energy and V = a 3 is the volume of the universe. Another important 
quantity is the entropy which is given by, 


S = 


_A_ 
4 G' 


(5.2) 
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Figure 4. Typical behavior of <5 in terms of t for finite A 2 (extended Chaplygin gas), n = 0 (blue), 
n = 1 (green), n = 2 (yellow), n = 3 (red). 


Using 8itG = 1 and A = An/H 2 [64] gives us the following expression for entropy in terms 
of the Hubble expansion parameter, 

S = 8t t 2 H~ 2 . (5.3) 


Using the relations (2.4) and (2.8) we have Hubble expansion parameter in terms of V, 
which is decreasing function and yields to a constant as expected, therefore we can extract 
the entropy in terms of V as follows, 


24vr 2 (y 10A2 e 37r - l) 

1 + V-l + 5U 10A 2e 3,r + V 10A2 e 3n ' 


(5.4) 


It is clear that the entropy is an increasing function of V (and therefore is an increasing 
function of time) which tells that the generalized second law of thermodynamics is valid 
as expected. These are obtained by choosing A 2 > 0 and we usually fix A 2 =4/3 to have 
Hi = 1/3 and H ~ 3 in agreement with some previous works such as [65]. 

Then, using the first law and the second law of thermodynamics, the temperature can be 
written as [60], 


(p + p)v 

S 


(5.5) 


Assuming Vo ~ 0.5 as initial volume of the universe gives the behavior of the temperature 
as illustrated in the Fig. 5. We can see high temperature at the early universe which 
decreases at the late time. It means that the temperature T> 0 as V —>• Vo and T —>• 0 as 
V Vo, therefore the third law of thermodynamics is satisfied for the extended Chaplygin 
gas model at the second order. 

It is interesting to obtain a relationship between the temperature and the energy density. 
Using the tan” 1 ( pecg + 1) ~ vr/2 approximation and above relations we can obtain, 


T = 


— 3-7T 

A 2 e 10 a 2 
247r 2 


(P 3 + P 2 


- 2 ) 


8 - 


p 2 + 2 p + 2 
p 2 - 2p + 2 


1 

10A 2 


(5.6) 
















while the exact relation yields to the behavior represented in the Fig. 6. Again we can 
see a good agreement between approximate solution and exact solution. As we expected 
the energy density increases by increasing temperature, and we find p oc T 1 / 3 at the early 
universe. 



Figure 5. The temperature in terms of V for Ai = 1.3. 



Figure 6. The energy density in terms of the temperature for A 2 = 1.3. Green line corresponds 
to the exact solution and blue line corresponds to the approximation solution. 


In order to verify the thermodynamic stability conditions of the extended Chaplygin 
gas along its evolution, it is necessary to have the following relations, 


dp 

dV 

p dS_ 

' dT 


< 0 , 


> 0 . 


(5.7) 
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The last equation is indeed the heat capacity which can be rewritten in terms of V [58], 

(5.8) 


d T 


C = V^\ — 
dV \dV 


Figs. 7 and 8 show that both conditions in (5.7) are satisfied without any restriction. It 
means that the extended Chaplygin gas behaves as a thermodynamically stable fluid (at 
least to the second order). 



v 


Figure 7. ^ in terms of V for A 2 = 1.3. 



V 


Figure 8. Heat capacity in terms of V for A 2 = 1.3. 

In order to get the partition function, we need to calculate free energy which is given 

by, 

F = U-TS , (5.9) 

where internal energy U is obtained using the equations (2.8) and (5.1). There is also 
another interesting quantity so called enthalpy (usually denoted by H but we denote with 
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77 to avoid conflict with Hubble parameter) which is given by, 


77 = U + pV. 


(5.10) 


In the plots of the Fig. 9 we represented the behavior of free energy and enthalpy in 
terms of V. We can see that the free energy is an increasing function of V, while enthalpy 
decreased to reach a constant at the late times. The enthalpy behaves as Hubble parameter 
with larger value at the initial time and smaller value at the late times. There is a critical 
V (close to the present epoch) where 77 = H. 



v 



Figure 9. Free energy (left) and enthalpy (right) in terms of V for A 2 = 1.3. 

Now, partition function can obtained using the following relation, 

Z = e" t. (5.11) 

In the Fig. 10, we can see behavior of the partition function with V. The maximum 
value of Z is near the early universe and it has small value at the present epoch. Before 
the extremum point the value of the free energy is negative, and the extremum point 
corresponds to zero free energy. 


6 Cosmological parameters 


We can investigate some cosmological parameters by varying thermodynamical quantities. 
An important parameter is the deceleration parameter given by, 


q = -1 + 


p + p 
2 H 2 ' 


( 6 . 1 ) 


Therefore we can study the behavior of the deceleration parameter in terms of V and T 
(see Fig. 11 ). As expected (7 —>- — 1 is obtained at the late times. Also, deceleration to 
acceleration phase transition happens in this model. We give a plot of the deceleration 
parameter in terms of redshift by Fig. 12 to see at what z the late acceleration begins for 
various possible values of A 2 . For instance, in the case of A 2 = 1.3 we see that universe 
starts accelerating at z < 0.09. 
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Figure 11. Deceleration parameter in terms of V (left) and temperature (right) 



Figure 12. Deceleration parameter in terms of redshift for A 2 = 0.3 (blue), An — 0.5 (cyan), 
An = 0.9 (green), A 2 = 1 (black) and A 2 = 1.3 (red). 
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The behaviour of the equation of state parameter u> = p/p is also illustrated in the 
Fig. 13. As expected it yields to -1 at the late times. Some time at the early universe we 
have oj = 0 so the extended Chaplygin gas behaves as dust (p = 0) and exchanges to de 
sitter at the late times. 



Figure 13. w in terms of V for A 2 = 1.3. 


Finally, we can investigate the squared sound speed (as promised before) which is an 
important parameter to study the stability of the model. As we know the fluid will be 
stable if, 


v l = 


dp 

= TT- > 0. 


dp 


( 6 . 2 ) 


In Fig. 14 we choose A 2 = 1.3 and see that C/ > 0 without any restriction, therefore there 
is no phase transition and the model will be stable at all times. 

So we demonstrated the stability of the ECG model using two separate ways, analyzing 
our model cosmologically and thermodynamically. Agreement of cosmological parameters 
with observational data supports our model. 


7 Conclusion 

In this paper, we considered extended Chaplygin gas at the second order which recover 
quadratic barotropic equation of state. First of all we obtained energy density in terms 
of scale factor and found an analytical expression for the scale factor. Then, we studied 
the evolution of density perturbations in both relativistic and Newtonian regime. In the 
relativistic regime we have shown that the perturbations may vanish for higher wavelength 
of perturbations. This is completely opposite with the modified Chaplygin gas where 
perturbations may vanish at small wavelength of perturbations. At the Newtonian limit 
we have shown that the perturbations vanish without any limitation. We also studied 
thermodynamical quantities of the model such as temperature and entropy. We confirmed 
that the thermodynamics laws satisfied and fluid is thermodynamically stable all times. 
Stability of the fluid also is investigated by using the sound speed. We found that there is 
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Figure 14. squared sound speed in terms of V for A 2 = 1.3. 

no phase transition and critical point in this model. Partition function is also calculated in 
the model using the free energy. Using the numerical analysis we found that the free energy 
is an increasing function of time (also V). Analyzing of enthalpy yields to an interesting 
result that is similar to the Hubble parameter at about present epoch. We investigated 
the behavior of some important cosmological parameters in terms of thermodynamical 
quantities. For example q —>• — 1 and u —l — 1 was shown to have desired behaviours at 
late times. We can also consider higher order terms and expect to obtain similar results. 
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